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Abstract 

We investigate transmission protocols for relay-assisted free-space optical (FSO) systems, when multiple 
parallel relays are employed and there is no direct link between the source and the destination. As alternatives to 
all-active FSO relaying, where all the available relays transmit concurrently, we propose schemes that select only 
a single relay to participate in the communication between the source and the destination in each transmission 
slot. This selection is based on the channel state information (CSI) obtained either from all or from some of the 
FSO links. Thus, the need for synchronizing the relays' transmissions is avoided and the slowly varying nature 
of the atmospheric channel is exploited. For both relay selection and all-active relaying, novel closed-form 
expressions for their outage performance are derived, assuming the versatile Gamma-Gamma channel model. 
Furthermore, based on the derived analytical results, the problem of allocating the optical power resources to the 
FSO links is addressed, and optimum and suboptimum solutions are proposed. Numerical results are provided for 
equal and non-equal length FSO links, which illustrate the outage behavior of the considered relaying protocols 
and demonstrate the significant performance gains offered by the proposed power allocation schemes. 
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I. Introduction 

The constant need for higher data rates in support of high-speed applications has led to the development 
of the Free Space Optical (FSO) communication technology. Operating at unlicensed optical frequencies, FSO 
systems offer the potential of broadband capacity at low cost 12, and therefore, they present an attractive 
remedy for the "last-mile" problem. However, despite their major advantages, the widespread deployment of 
FSO systems is hampered by major impairments, which have their origin in the propagation of optical signals 
through the atmosphere. Rain, fog, and atmospheric turbulence are some of the major atmospheric phenomena 
that cause attenuation and rapid fluctuations in the received optical power in FSO systems, thereby increasing 
the error rate and severely degrading the overall performance (U. 

In the past, several techniques have been applied in FSO systems for mitigating the degrading effects of 
the atmospheric channel, including error control coding in conjunction with interleaving 0, multiple-symbol 
detection (4]j, and spatial diversity 0-1171. Among these techniques, spatial diversity, which is realized by 
deploying multiple transmit and/or receive apertures, has been particularly attractive, since it offers significant 
performance gains by introducing additional degrees of freedom in the spatial dimension. Thus, numerous FSO 
systems with multiple co-located transmit and/or receive apertures, referred as Multiple-Input Multiple-Output 
(MIMO) FSO systems, have been proposed in the technical literature J5]-||7]. However, in practice, MIMO 
FSO systems may not always be able to offer the gains promised by theory. This happens in cases where the 
assumption that all the links of the MIMO FSO system are affected by independent channel fading becomes 
invalid |7]. Furthermore, since both the path loss and the fading statistics of the channel are distance-dependent, 
a large number of transmit and/or receive apertures is required in long-range links in order to achieve the desired 
performance gains, thus increasing the complexity of MIMO FSO systems. 

In order to overcome such limitations, relay-assisted communication has been recently introduced in FSO 
systems as an alternative approach to achieve spatial diversity ||8l- lfTTll . The main idea lies in the fact that, by 
employing multiple relay nodes with line-of-sight (LOS) to both the source and the destination, a virtual multiple- 
aperture system is created, often referred as cooperative diversity system, even if there is no LOS between the 
source and the destination. In JSJ, various relaying configurations (cooperative diversity and multihop) have 
been investigated under the assumption of a lognormal channel model. Subsequently, several coding schemes 
for 3-way cooperative diversity FSO systems with a single relay and a direct link between the source and the 
destination have been proposed in J9), while the performance of such systems has been investigated in ifTOll and 
ifTTTl assuming amplify-and-forward and decode-and-forward relaying strategies, respectively. It is emphasized 
that in all these previous works, all the available relays participated in the communication between the source 
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and the destination, requiring perfect synchronization between the relays such that the FSO signals can arrive 
simultaneously at the destination, while the optical power resources are equally divided between all FSO links. 

In view of the above, in this paper we present alternative transmission protocols which can be applied to 
relay-assisted FSO systems with no LOS between the source and the destination. For the signaling rates of 
interest, the atmospheric channel does not vary within one packet. Thus, channel state information (CSI) can 
be easily obtained for all or for some of the involved links. Capitalizing on this fact, the presented protocols 
select only a single relay to take part in the communication in every transmission slot, thus avoiding the need 
for synchronization between the relays. It should be noted that similar relay selection protocols have been 
also proposed in the context of radio-frequency relaying systems |[l2l - |[T4l . In particular, two types of relay 
selection protocols are presented: the select-max protocol that selects the relay that maximizes an appropriately 
defined metric and requires CSI from all the available FSO links, and the distributed switch and stay (DSSC) 
protocol which switches between two relays and requires CSI only from the FSO links used in the previous 
transmission slot. Furthermore, assuming the versatile Gamma-Gamma channel model O and decode-and- 
forward relay nodes, we derive novel closed-form analytical expressions for the outage performance of the 
proposed transmission schemes, as well as the scheme where all the available relays transmit simultaneously; 
thus, extending the analysis presented in [8 ] to the case of the Gamma-Gamma channel model. Finally, based on 
the derived outage results, we address the problem of optimizing the allocation of the optical power resources 
to the FSO links for minimization of the probability of outage; hence, rendering the relay-assisted FSO system 
under consideration more power efficient. 

The remainder of the paper is organized as follows. In Sectionim the system model and the considered relaying 
protocols are discussed. The outage performance of the relaying protocols under investigation is analyzed in 
Section |III1 while the problem of optimizing the allocation of the optical power resources to the FSO links is 
addressed in Section JV] Numerical results for various relay-assisted FSO architectures are presented in Section 
|VJ and, finally, concluding remarks are provided in Section |VlJ 

II. System Model 

The system model under consideration is depicted in Fig. [T] In particular, we consider an intensity-modulation 
direct detection (IM/DD) FSO system without LOS between the source, S, and the destination, D, and the 
communication between these two terminals is achieved with the aid of multiple relays, denoted by Ri, 
i 6 {1, ...,N}. The source node is equipped with a multiple-aperture transmitter, with each of the apertures 
pointing in the direction of the corresponding relay, and an optical switchQ, which either allows the simultaneous 

'Optical switches can be implemented with either spatial light modulators (SLM) 1151 Ch. (27)] or optical MEMS devices 1161 . 
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transmission from all the transmit apertures or selects the direction of transmission by switching between the 
transmit apertures. 

The presence of a large field-of-view (FOV) detector at the destination is assumed allowing for the simultane- 
ous detection of the optical signals transmitted from each relay. Moreover, all optical transmitters are equipped 
with optical amplifiers that adjust the optical power transmitted in each link. The relaying terminals use a 
threshold-based decode-and-forward (DF) protocol; that is, they fully decode the received signal and retransmit 
it to the destination only if the signal-to-noise ratio (SNR) of the receiving FSO link exceeds a given decoding 
threshold. Finally, throughout this paper, we assume that binary pulse position modulation (BPPM) is employed. 

A. Signal and Channel Model 

For an FSO link connecting two terminals A and B, the received optical signal at the photodetector of B is 
given by 

r s nTh ioAnP+h . „ + Ph) + n s 

(1) 



?B 



r]T b (p AB P t h AB + P b ) + n s 
r]T b P b + n n 

where r s and r n represent the signal and the non-signal slots of the BPPM symbol, respectively, while pAsPt 
and P b denote the average optical signal power transmitted from A and the background radiation incident on 
the photodetector of B, respectively. Furthermore, pab represents the percentage of the total optical power P t 
allocated to the FSO link between terminals A and B, h AB is the channel gain of the link, r\ is the photodetector's 
responsivity, T b is the duration of the signal and non-signal slots, and n s and n n are the additive noise samples 
in the signal and non-signal slots, respectively. Since background-noise limited receivers are assumed, where 
background noise is dominant compared to other noise components (such as thermal, signal dependent, and 
dark noise) 0, ifTTl , the noise terms can be modeled as additive white Gaussian, with zero mean and variance 
a n = ^~2~- After removing the constant bias r[T b P b from both slots, the instantaneous SNR of the link can be 
defined as H 

7 - = N • (2) 

Due to atmospheric effects, the channel gain of the FSO link under consideration can be modeled as 

h AB=h AB h AB (3) 

where Kab accounts for path loss due to weather effects and geometric spread loss and h AB represents irradiance 
fluctuations caused by atmospheric turbulence. Both h AB and h AB are time-variant, yet at very different time 
scales. The path loss coefficient varies on the order of hours while turbulence induced fading varies on the order 
of 1-100 ms HI. Thus, taking into consideration the signaling rates of interest, which range from hundreds 
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to thousands of Mbps, the channel gain can be considered as constant over a given transmission slot, which 
consists of hundreds of thousands (or even millions) of consecutive symbols. 

The path loss coefficient can be calculated by combining the Beer Lambert's law [2 1 with the geometric loss 
formula CD pp. 44], yielding 



h. 



■ exp (—vdj. 



(4) 



' (D T + e T d AB ) 2 

where Dr and Dt are the receiver and transmitter aperture diameters, respectively; 9t is the optical beam's 
divergence angle (in mrad), d AB is the link's distance (in km), and v is the weather dependent attenuation 
coefficient (in 1/km). 

Under a wide range of atmospheric conditions, turbulence induced fading can be statistically characterized 
by the well-known Gamma-Gamma distribution [2]. The probability density function (pdf) for this model is 
given by 

a AB + ,) i " 

2K B D ; 



-x 



(5) 



T(a AB )T((3 Al 

where T (•) is the Gamma function lfT8l Eq. (8.310)] and K u (•) is the uth order modified Bessel function of the 
second kind lfT8l Eq. (8.432/9)]. Furthermore, aAB and (5ab are parameters related to the effective atmospheric 



conditions via a t 



cxp 



0A9a 2 R 



and P, 



exp 



0.51oj 



l+0.69o- R 5 



12 \T 



-1 



0, where 



a\ = 1.23C 2 (x) 6 ^Ib denotes the Rytov variance^, C% is the weather dependent index of refraction structure 
parameter, and A represents the wavelength of the optical carrier. 



B. Mode of Operation 

Throughout this work, three different cooperative relaying protocols are considered: the all-active protocol, 
originally presented in ||8], where all the available relays are activated, and the select-max and the distributed 
switch and stay combining (DSSC) protocols, which are both based on the concept of selecting a single relay. 

1 ) All-active: In this relaying scheme, the source activates all relays and the total power is divided between 
all available FSO links. Since the relay nodes operate in the DF mode only the relays that successfully decode 
the received optical signal remodulate the intensity of the optical carrier and forward the information to the 
destination. At the destination, owing to the presence of a large FOV aperture, aperture averaging occurs |]6] and 

2 The Rytov variance is indicative of the strength of turbulence-induced fading. More specifically, values tj\ < 1 correspond to weak 
turbulence conditions, while values o\ > 1 correspond to the moderate-strong turbulence regime [1191 . 
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all the received optical signals are added. Hence, assuming perfect synchronization, the output of the combiner 
can be expressed as 



VmgD 



(6) 



V T b P b + n n 

where D denotes the decoding set formed by the relays that have succesfully decoded the signal. Since the 
total power is divided between all available links, it follows that £ i=1 {Psr + Pr o) = !• 

The advantage of this scheme is that CSI is not required neither at the transmitter nor the receiver side, since 
the source transmits to all available relays, regardless of their channel gain. However, since it is assumed that 
all the signals arrive at the destination at the same time, this scheme requires accurate timing synchronization 
in order to account for the different propagation delays of the different paths, resulting in high complexity. 

2) Select-Max: This relaying protocol selects a single relay out of the set of N available relays in each 
transmission slot. In particular, the relay which maximizes an appropriately defined metric is selected. This 
metric accounts for both the S-R4 and R t -D links and reflects the quality of the zth end-to-end path. Here, we 
adopt the minimum value of the intermediate link SNRs, 

7i = min(7 SR .,7 R . D ) , (7) 

as the quality measure of the zfh end-to-end path, which will be referred as the "min equivalent SNR" throughout 
the paper. Note that $7} represents an outage-based definition of the selection metric, in the sense that an outage 
on the ith end-to-end link occurs if 7$ falls below the outage threshold SNR. Hence, the single relay that is 
activated in the select-max relaying protocol, R b , is selected according to the rule 

b = argmax 7$ . (8) 
te{l,..j\r} 

Since a single relay is activated in the select-max protocol, the total available optical power is divided between 
the S-R b and R b -D links, i.e., p SR + p R D = l,and in the case that R b has successfully decoded the received 
optical signal, i.e., b € D, the signal at the destination can be expressed as 

riT b (p RbD P t h RbD +P b )+n s 
r]T b P b + n n 

This relaying scheme requires the CSI of all the available S-Ri and Ri-D FSO links in order to perform the 
selection process. This can be achieved by some signalling process that takes advantage of the slowly-varying 
nature of the FSO channel . Here, each receiver estimates the correspponding link CSI and feeds it back to the 
source through a reliable low-rate RF feedback link. It is emphasized that since only one end-to-end path is 
activated in each transmission slot, only one signal arrives at the destination and thus synchronization between 
the relays is not needed. 



(9) 
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3) DSSC: Requiring less CSI than select-max, the DSSC protocol applies to the case where there are only 
two relays available and one of them is selected to take part in the communication between the source and 
the destination, in a switch-and-stay fashion [13]. More specifically, in each transmission slot the destination 
compares the min equivalent SNR of the active end-to-end path with a switching threshold, denoted by T. 
If this SNR is smaller than T, the destination notifies the source and the other available relay is selected for 
taking part in the communication, regardless of its end-to-end performance metric. 

Mathematically speaking, denoting the two available relays by R\ and R 2 and the min equivalent SNR of 
the ith end-to-end path during the jth transmission period by 7? , the active relay in the jth transmission period, 
R 3 b , is determined as follows: 

. , . [ Ri when 4 > T 

if R{~ 1 = Ri then R{= { (10) 

! R 2 when j{ < T 

and 

, I R2 when ~i'i>T 
if R{~ 1 = R 2 then R{= { . (11) 

[ R x when 7^ < T 

Hence, in the case that R 3 b has successfully decoded the received signal, the optical signal at the destination is 
given by 

rffi, ( p , P t h , + P h ) + n s 

(12) 



rjT b P b + n n 

Since in this protocol only a single relay assists in the communication between the source and the destination, 
the power allocation rule of the select-max protocol also holds for DSSC relaying. 

When there are more than two available relays in the system, i.e., N > 2, a modified version of DSSC 
protocol could initially sort all the available paths based on their end-to-end distance, defined as 

di = max(d SR .,d R . D ) (13) 

with i = l,...,N, and, then, use as R\ and R 2 the two relays that correspond to the paths with the minimum end- 
to-end distance. It should be noted that end-to-end distance is an indicative of the path's end-to-end performance, 
taking into consideration that both path loss and Rytov variance are monotonically increasing with respect to 
the link distance. 

The simplicity of this scheme compared to the select-max protocol lies in the fact that only the CSI of 
the active end-to-end path is required for the selection process, resulting in less implementation complexity. 
Furthermore, as in the select-max scheme, no synchronization among the relays is needed, since only one 
end-to-end path is activated in each transmission slot. 
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III. Outage Analysis 
At a given transmission rate, Vq, the outage probability is denned as 

Pout(r ) = Pr{C( 1 ) <r }, (14) 

where C (■) is the instantaneous capacity, which is a function of the instantaneous SNR. Since C (•) is 
monotonically increasing with respect to 7, (fl4l) can be equivalently rewritten as 

iWr )=Pr{7<7th}, (15) 

where 7^ = C _1 (7*0) denotes the threshold SNR. If the SNR, 7, drops below 7^, an outage occurs, implying 
that the signal cannot be decoded with arbitrarily low error probability at the receiver. Henceforth, it is assumed 
that the threshold SNR, 7^, is identical for all links of the relaying system. 

A. Outage Probability of the Intermediate Links 

Since DF relaying is considered, an outage event in any of the intermediate links may lead to an outage of 
the overall relaying scheme. Therefore, the calculation of the outage probability of each intermediate link is 
considered as a building block for the outage probability of the relaying schemes under investigation. 

By combining Q with ([TBI ), the outage probability of the FSO link between nodes A and B is defined as 

(r] 2 T 2 p 2 R P 2 h 2 1 
Pout,AB = Prj 1 hPA ^ 1 < lth\ (16) 

which can be equivalently rewritten as 

Pout,AB = Pr [h AB < = (17) 

I Ii ab pabPm J 

where Pu is the power margin given by Pm = l]f Pt ■ Using the cumulative density function (cdf) of the 
Gamma-Gamma distribution EDI Eq. (7)], the outage probability of the FSO link between nodes A and B can 
be analytically evaluated for any chab and (3ab, yielding 

1 

where [•] is the Meijer's G-function US Eq. (9.301)]. 

To gain more physical insights from (fLSl l. it is meaningful to explore the outage probability in the high power 
margin regime. 

Theorem 1: For high values of power margin and when (a AB — f3 AB ) $ Z, the outage probability of the FSO 
link between nodes A and B can be approximated by 

( " A nP A n \ q A» 

d ^ l{j> ab — q AB ) \h AB P M pABj . 1Q . 
rout,AB — tt? v r (a — r K^) 

AB/ \r AB / QAB 
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where p AB = max (a AB , f3 AB ) and q AB = min (a AB ,/3 AB ). 

Proof: A detailed proof is provided in Appendix I. ■ 
It should be noted that in the analysis that follows it is assumed that (a AB — f3 AB ) ^ Z holds for every 
possible FSO link. Although this condition may seem restrictive, it can be relaxed in practical applications by 
inserting an infinitely small perturbation term e, so that (a AB — j3 AB + e) ^ Z, when (a AB — f3 AB ) £ Z [21]. 



B. Outage Probability of All-Active Relaying 

In this scheme an outage occurs when either the decoding set D is empty or the SNR of the multiple-input 
single-output link between the decoding relays and the destination falls below the outage threshold. Hence, the 
outage probability of this scheme can be expressed as [H] Eq. (30)] 

2 ( 1 1 

E PB mB h amB <-p-)Pr{S(n)} 



out 



71=1 



(20) 

I'M II 

where S (n) denotes the nth possible decoding set, 2^ is the total number of decoding sets, and Pr {S (n)} is 
the probability of event {D = S (n)} given by 



Pr{S(n)} 



n pv bs Rm >ith} n pr {7 SRm <7^} 

1 



n 



l-PrU sflm < p 

P SR m M 



n pr {^ 



< 



i 



P S R m P M 



(21) 



1 ) Exact Analysis: In order to evaluate (l20l . the cdf of the sum of weighted non-identical Gamma-Gamma 
variates, h s , n) = Yl P h RmI} , needs to be derived first. However, to the best of the authors' knowledge, 

m£S(n) m 

there are no closed-form analytical expressions for the exact distribution of the sum of non-identical Gamma- 
Gamma variates. Therefore, the numerical method of ||22] Eq. (9.186)], which is based on the moment generating 
function (MGF) approach, is applied and thus the cdf of h g[n) , denoted as Fh ai (•), is evaluated via 



F h fx) 



2-^exp(4) ^ 




2 



r n {M BmD 

mSS(n) 



V 



4)) 



A_ 

2x 



L+k 

+ £(-l)'Re 

1=1 



n (m 

mSS(n) V 



A+j2nl 
2x 



A+j2irl 
2x 




(22) 



where Ai, 



is the MGF of the channel gain of the R m D FSO link given by [23] Eq. (4)], while the 



parameters A, K, L are calculated based on the numerical error term obtained by [22, Eq. (9.187)]. 
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Theorem 2: The outage probability of the all-active relaying protocol in Gamma-Gamma fading is given by 

9™ / 



P 



out 



e n 

n=l ragS(n) 



1 



1 



V 



r(a SRm )r(p 3Rm ) 



G 



2.1 

1,3 



SR m ^SR„ 



Ps Rm P M 



KJr(/3 SR J 



G 



2,1 
1,3 



S(») I R 



(23) 



Proof: The proof follows straightforwardly by combining (120} with (TT8T ) and (l22l . ■ 
2) Asymptotic Analysis: In order to gain more physical insights into the performance of the relaying protocol 

under consideration, we further consider the high power margin regime, i.e., when Pm — > oo. In order to perform 

this analysis, an asymptotic expression for Fh (•) needs to be derived first. 

Lemma 3: For high values of power margin, the cdf for the weighted sum of non-identical Gamma-Gamma 

variates that corresponds to decoding set S (n), h s(n) , can be approximated as 



F h (x) 



n 

meS(n) 



a R m ol 3 R 7 , 
^r dPr„ 



r (lR m nF(p R „ 



D "R m D , 



r K mD ) r (^. mD ) 



E 9 RmD 



(24) 



E 

im£S(n) 



E q Rm n 

v mS5(n) 



Proof: A detailed proof is provided in Appendix II. ■ 
The asymptotic expression for the outage probability of the all-active relaying scheme is given by the following 
theorem. 

Theorem 4: For high values of power margin, the outage probability of the all-active relaying scheme in 
Gamma-Gamma fading can be approximated by 

I I SRm^SRm 1 \ 



P 



out 



e n 

n=l m<£S{n) 



n 

mSS(n) 



V 



17 (" S « j r (/^K j VSR 



7 



R m opR m n\ qRmD r(q RmD )r(p RmD -q RmD ) 



h R m nP Rm o 



E Q RmD r E ?«, 



r K mD ) r (^ mD ) f i 
p^ 



E Q S R m + E <? RmD 



• (25) 



meS(n) 



mS5(n) 



Proof: We first observe that in the high power margin regime, i.e., Pm — ► oo, (TUT) can be approximated 



by 



Pr{S»}« [J Pr(/ isRm < 1 —\ 

7a, \ I Ps R PM J 



(26) 



m^S{n) 
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Hence, by combining (1201 ) with ( |T9l ), (1241 ). and (1261 ). the asymptotic expression in (1251 ) is obtained. This concludes 
the proof. ■ 

An important result derived from the asymptotic expression of the previous theorem is the diversity gain of 
the transmission protocol under consideration which is summarized in the ensuing corollary. 

Corollary 1: For the all-active relaying protocol, the diversity gain of a relay-assisted FSO system with N 
relays is given as 



Gd 



mm 

n£{l 2"} 



,m^S(ra) meS(n) 



Proof: We define the diversity order as 

G d 



lim 



log Pout 



(27) 



(28) 



Pm^oo log Pm 

Hence, (l27l) follows straightforwardly by observing that the term which corresponds to the power of Pm with 
the minimum of ^ Qsn m + Yl 9n m r> dominates in the sum of (|25T ). when Pm — > oo. ■ 

\m<tS(n) m meS(n) / 

C. Outage Probability of Select-Max Relaying 

In the select-max protocol a single relay out of the N available relays is selected according to the selection 
rule in ©. Hence, the outage probability of the relaying scheme under consideration is given by 



N 



Pom = Pom {Ri n ... n R N } = Yl Pout {Rb} 



(29) 



6=1 



where P out {Rb} denotes the probability of outage when only relay R^ is active. Given that R}, is active, an 
outage occurs when either Rt or D have not decoded the information successfully, i.e., 



Pout {Rb} 



Pr {(7s« 6 <lth)u( lRbD <jth)} 



l 



1 - 1 - Pr ih gK < - 

V " h sn h Psn h 



(30) 



Hence, the probability of outage of the select-max relaying scheme is obtained by combining d29l with d30l ). 

1 ) Exact Analysis: The following theorem provides an accurate analytical expression for the performance 
evaluation of the select-max relaying scheme. 

Theorem 5: The probability of outage of a relay-assisted FSO system that employs the select-max relaying 
protocol in Gamma-Gamma turbulence-induced fading is given by 



N 



out 



n > 



i 



i 



6=1 



X 1 



r K) r (^ 



r«2,i 

"^1,3 



hsR.PsnPM 



a SR h ' fisR h ' 



2,1 



r(a„ D )r(/3„ D ) Gl > 3 



0l R i D^R b D 

h Rh nP Rh n P M 



^RhD 1 Pr,,D 1 



(31) 
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Proof: The proof follows straightforwardly by combining (129) with (1301 ) and (fT8l ). ■ 
2) Asymptotic Analysis: In order to gain more physical insights into the performance of the relaying protocol 
under consideration, we investigate its asymptotic behavior when Pm — > oo, in the ensuing theorem and 
corollary. 

Theorem 6: For high values of power margin, the outage probability of the select-max relaying scheme can 
be approximated as 

/ r (PsR b -1 S R b ) ( a SR b P S R b \ q sn„ T {PR b n-lR b n) / a R b pP Rb n \ 



P, 



out 



N 

n 

6=1 



V 



r(a SRb )r(/3 s 



PsrPm 



+ 



(32) 



Proof: The proof starts by observing that as Pm —> oo the probability of outage given relay Rt, is active, 
can be approximated by 



Pout {Rb} « Pr ^ h SRh < 



+ Pr\h RbD < - Pm 

n R h nPR h n 



(33) 



Hence, by combining d29l ) with ( fT9l and (l33l . the asymptotic expression in (l32l is obtained. This concludes 
the proof. ■ 
Corollary 2: The diversity gain of a relay-assisted FSO system employing the select-max relaying protocol 
and N relays is given as 

N 

G d = Y^ min (q S R b >q RbD )- (34) 

6=1 

Proof: The proof follows straightforwardly from (|32~V When P^/ — > oo> the term that corresponds to the 
power of Pm with the minimum of (q 3R ,q R D ) dominates in the sum inside the product. Hence, after taking 
the product of the dominating terms, the diversity order is obtained. ■ 

D. Outage Probability of DSSC Relaying 

In the DSSC protocol, the selection of the single relay which takes part in the communication is based on 
(TTOb and (fTTl ). Hence, an outage occurs when there is an outage either in the end-to-end link of the first relay, 
given that the first relay is selected in the jth transmission slot, or in the end-to-end link of the second relay, 
given that the second relay is selected in the jth transmission slot, i.e., 

P out = Pr { [W b = Px) n (tJ < 7th) } + Pr { (p j 6 = P 2 ) fl (4 < lt h) } (35) 

Following the analysis of ll22l Sec. (9.9.1.2)], the above equation can be rewritten as 

' ffffff + ^y<j? } (Pr {71 < 7th} + Pr { 7 2 < 7*4) , 7th < T 

It^^Si ( p r {71 < 7th} + Pr {72 < 7th} - 2) (36) 



P. 



out 



+ 



(Pr{ 7l < 7th } Pr{ 72 <T}+Pr{ 7l <T} Pr{ 72 <7 th }) 
Pr{ 7 i<T}+Pr{ 72 <T} 



, 7th > T 
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1 ) Exact Analysis: The following theorem provides an accurate analytical expression for the performance of 
the DSSC relaying scheme. 

Theorem 7: The probability of outage of a relay-assisted FSO system that employs the DSSC relaying 
protocol is given by 

T<P M 



P 



out 



[ MiK(|) ( Fh (Jj\+ Fh (+- 

+ F hl (i) + F h3 (i) • T>P ™ 



(37) 



where 



1-1 



x 1 



r(« Sfl Jr(/3 SB J 



r«2,i 

^1,3 



-G 



2,1 
1,3 



^ RiD 1^ R i D 



1 
1 

^Hi D ! fin. D 1 



(38) 



and f = ?m 



N T 



Proof: We first note that the cdf of the min equivalent SNR defined in |7]) can be expressed according to 
(21 pp. 141] as 

F 7, (*) = 1 - (l - F lsRi (x)) (l - F lRiD (x)) (39) 
which is equivalent, after a variable transformation, to 

T. \ \ i 

1 — Pout,RiD 



F 7< (x) = 1 - 1 - P c 



out,SRi 



P Ri oVTbPt 



(40) 



.PsR^TbPt, 

Using (fT6l ) the cdf of the min equivalent channel gain is derived as (|38l ) and, thus, according to d36l ), d37l ) is 
obtained. This concludes the proof. ■ 
Corollary 3: The performance of the DSSC relaying scheme is minimized when T = Pm and in that case 
it becomes equal to that of the select-max scheme for two relays, R\ and i?2- 

Proof: Following the analysis in ll22l (Ch. 9.9.1.1)], the performance of DSSC relaying is minimized when 
T = Pm is set. In that case, d37l) yields 



Pout = 11 (1 " (1 " Pout, SRi ) (1 " Pout^D)) 



(41) 



4 = 1 



which is equivalent to (f3TT > when N = 2. This concludes the proof. 
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2) Asymptotic Analysis: In order to gain more physical insights of the DSSC protocol with optimized 
threshold, we investigate the asymptotic behavior of its performance when Pm — > oo in the ensuing corrolary. 

Corollary 4: The minimum outage probability of the DSSC relaying protocol can be approximated at the 
high power margin regime, by (l32i with N = 2, and the maximum achieved diversity gain is given by 

2 

G d = Y^ min (<l SRi , tfiv, ) ■ ( 42 ) 
i=i 

Proof: The proof follows straightforwardly by combining (|4H with (l32l ) and (l34l ). ■ 

IV. Optimal Power Allocation 

In this section, we are interested in optimizing the optical power resources in both the S-Ri and R{-D links, 
in order to minimize the outage probability of the relay-assisted FSO system for a given total optical power. 
Hence, in the following, we optimize the parameters p SR and p R D for each of the relaying protocols under 
consideration. 



A. Power Allocation in All-Active Protocol 

Since in the all-active scheme the power is divided among all the underlying links, the minimization of its 
outage probability is subject to two constraints; the total power budget of all links is equal to P t and the optical 
power emitted from each transmitter is less than P t . Consequently, the optimum power allocation can be found 
by solving the following optimization problem 

min P out 

Em=l {pSR m + /?r to d) = 1 ^) 

subject to < o < p SRm < 1, m = 1, ..JV 
^ 0<p RmD <1, m = l,...N 
where P out is given by (l23l ). It should be noted that that the above optimization problem is convex problem. 
This can be explained as follows. Since the objective function is an outage probability, it is convex according 
to ||25l . Furthermore, since all the constraints are linear, they form a convex set 11261 , which leads to a convex 
optimization problem and, thus, a unique optimal solution. 

Using the exact outage expression in (|23l ), it is difficult to find the optimum solution for the problem in (|43l i, 
even with numerical methods, due to the involvement of the Meijer's G-functions. Therefore, the asymptotic 
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expression of (125b is used as objective function instead and hence the optimization problem is reformulated as 

.1 ■ - j i / 

E </„ mD )r( E g, n ) ^S(n) 



mm 



2 N 

E 

71=1 



V 



x n 

mgS(n) 



h R m D p R m D 



(44) 



subject to < 



r K mD ) r (/ ? R mD ) 

0<p Sflm <1, m = l,...N 

^ 0<P flmD <1, m = l,..JV 
which is a geometric program that can be numerically solved using numerical optimization techniques, such as 

the interior point method [26 Sec. 14]. 

Since the derivation of the exact solution is cumbersome and motivated by the dependence of the outage 
probability on the link distance, the following suboptimal power allocation scheme for all-active relaying is 
proposed. 

Proposition 1: For all-active relaying, the fraction of the total optical power which is allocated to each link 
is given by 

P SR , = — and P R ,o = — (45) 



JV 

E (d SRm +d RmD ) 

m=l 

for the S-R{ and Ri-D links, respectively, with i = 1, N. 



E { d s Rm + d Rm n) 
m=l 



B. Power Allocation in Select-Max Protocol 

Similar to the all-active scheme, the outage probability of the select-max protocol can also be minimized by 
optimizing the optical power resources which are allocated to each of the links. However, in this scheme, both 
the objective function and the constraints are different from the problem in (|43T ). 

Since the total optical power is divided only between the S-Rb and Rb-D links of the active relay, the 
problem can be formulated as 

min P out {R b } 

[Ps Rb +P RhD =l (46) 
subject to I < p SRb < 1 

0<p R&D <l 

where P out {Rb} is the the probability of outage when Rb is active, given by (l30l ). Based on the same reasoning 
as in the previous relaying scheme, we conclude that the above optimization problem is convex, thus, leading 
to a unique optimal solution. 
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Due to the involvement of the Meijer's G-functions, it is again difficult to find the optimum solution if the 
exact expression in (l30l is used as objective function, even with numerical methods. Therefore, the asymptotic 
expression in (|33l ) is employed and hence the power allocation optimization problem is reformulated as 



r P 



mm 



SR h ISRu 



P SR -Pa. 



subject to < 



+ 



r (v) r K 



D )1R h D 



(47) 



< P SR , < 1 
0<p RbD <l 

Theorem 8: The power allocation parameters that minimize the outage probability of a relay-assisted FSO 
system when operating in Gamma-Gamma turbulence induced fading and employing select-max relaying, are 
given by 

i i 

(48) 



for the S-Rb and Rb-D links respectively, where b = 1,...,N, 5 SR 

r( P R bD -lR bD ) ( «R bD PR b ,, ' ''" ' ^ 



T {PsR b -1 S R b ) ( <*SR b Ps 



h R h D P ^- 



and to £ 



0, min 



^ S R b )r^ SRb 
is the unique real positive root of 



. 6 R 

R bD 



S(t) = 6 s s R * b+1 t q ™> +1 + 6 R n b b D D+1 t q » b ° +1 - 1. 



Proof: We first define the Langrangian associated with the optimization problem of (l47"T ) as 

^SJJ,, / J \ lR b D 



J = s s 



1 



+ 5 K 



,PsR b J \ P R b D 

Applying the method of Langrange multipliers, setting 



(49) 



(50) 



and 



0, and using the equality 



9 P S R b d pR b D 

constraint in (|47l , it is straightforward to show that the optimum power allocation coefficients are given by 
(l48l) . where to = — \ is the root of S (t), which has to take values in the interval 0, min ( j-^—, -r^—) in 

L \°SR b °R bD J\ 

order to satisfy the inequality constraints in (l47"V It should be noted that S (t) has a unique positive root which 
lies in this interval; this can be proved by applying the intemediate value theorem of continuous functions that 



0, min 



S SR ' S R 
bR b R bD 



and the Descartes' rule 



shows that S (t) has at least one real positive root in the interval 
of signs [27 ] which shows that this root is unique. ■ 
In order to avoid finding the root of the polynomial S (t) and motivated by the dependence of the outage 
probability on the link distance, the following suboptimal power allocation scheme for select-max relaying is 
proposed. 



January 21, 2013 



DRAFT 



SUBMITTED TO THE IEEE TRANSACTIONS ON COMMUNICATIONS 16 

Proposition 2: The fraction of the total optical power which is allocated to each link when the select-max 
relaying scheme is employed is chosen as 

p SRb = and p RbD = b — (51) 

for the S-Rb and Rb-D links, respectively, with b = 1, N. 



C. Power Allocation in DSSC Protocol 

Since a single relay is activated in each transmission slot by the DSSC protocol, the optimum power allocation 
scheme is obtained by minimizing the outage probability of the active end-to-end path. Hence, the optimization 
problem that has to be solved in this case is equivalent to the problem in (|47T ) and, therefore, the optimum and 
the suboptimum power allocation schemes of the select-max protocol can also be employed for DSSC relaying. 

V. Numerical Results and Discussion 

In this section, we illustrate numerical results for the outage performance of the considered relaying protocols, 
using the derived analytical expressions. In the following, we consider a relay-assisted FSO system with A = 
1550 nm and transmit and receive aperture diameters of Dr = Dt = 20 cm. Furthermore, we assume clear 
weather conditions with visibility of 10 km, which correspond to a weather-dependent attenuation coefficient 
of v = 0.1 ^ and an index of refraction structure parameter of C% = 2 x 10 -14 m~3. 

Fig. |2] depicts the outage performance of the presented relaying protocols for various numbers of relays, when 
the link distance is identical for all S-R4 and R{-D links and the optical power is equally divided between the 
active relays. Specifically, analytical results for the outage probability of a relay-assisted FSO system with a link 
distance of 2 km are plotted, as functions of the power margin for N = 2, 3, 4 relays, using the exact and the 
asymptotic outage expressions for each of the considered relaying protocols. We assumed p SR = p R D = ^ 
for the all-active and p SR . = p RO = ^ for the select-max and DSSC protocols, respectively. As benchmarks, 
Monte Carlo (MC) simulation results and the performance of an FSO system with N = 1, which is independent 
of the employed relaying protocol, are also illustrated in Fig. |2] As can be observed, there is an excellent match 
between simulation and analytical results for every value of N, verifying the presented theoretical analysis. 
Moreover, it is obvious that the select-max relaying scheme has a better performance compared to the all-active 
scheme in every case examined (performance gains of 2, 4, and 5 dB are observed for N = 2, 3, and 4, 
respectively). This result is intuitively pleasing, since the select-max protocol selects in each transmission slot 
the best end-to-end path out of the N available paths and allocates the total available optical power only to 
this path. Furthermore, when increasing the number of relays in the select-max and all-active protocols, it is 
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observed that the outage performance is significantly improved with respect to the single relay FSO system. In 
contrast, although the DSSC scheme with the optimum threshold offers significant performance improvement 
for N = 2 (its performance is identical with the select-max performance of N = 2), it remains unaffected by 
the increase of the number of relays. 

Fig. [3] depicts the outage performance of a relay-assisted FSO system employing the presented protocols 
and assuming different distances for each of the S-Ri and Ri-D FSO links. Specifically, two different system 
configurations are investigated. In the first system configuration, N = 2 and the link distances are given by 
vectors d SR = {2,1.5} and d RD = {1,2.5}, with the elements of the vectors respresenting the distances (in 
km) of the S-Ri and Ri-D links respectively, while in the second configuration N = 3, and the link distances 
are given by d SR = {2,1.5,1} and d RD = {1,2.5,3}. Fig. [3] reveals that the select-max relaying scheme 
offers significant performance gains compared to the all-active scheme, also for non-equal link distances. In 
particular, in the first configuration a gain of 2.5 dB compared to the all-active scheme is offered, while in 
the second configuration the offered gain is 3 dB. Furthermore, it can be easily observed that although in the 
second configuration the number of relays has been increased and the performance of both all-active and select- 
max relaying has been improved, DSSC with optimized threshold remains unaffected by this increase and its 
performance remains identical with the performance for the first configuration. This was expected, since DSSC 
uses only two end-to-end paths (those with the minimum end-to-end distance) and, therefore, the addition of 
extra paths with larger end-to-end distance will not improve the performance of this protocol. Finally, we note 
that simulation and analytical results are again in excellent agreement. 

Fig. @] illustrates the effect of power allocation in relay-assisted FSO systems employing the relaying proto- 
cols under consideration. Specifically, the performaces of the optimum and the proposed sub-optimum power 
allocation schemes, obtained by solving (l44b . (|47T) and from the empirical rules of (l45l . (1511 . respectively, are 
presented along with the equal power allocation, when the second system configuration of Fig. [3] is considered. 
It is obvious from Fig. @] that optimized power allocation offers significant performance gains compared to 
equal power allocation, irrespective of the employed relaying protocol. This was expected, since both the path 
loss and turbulence strength are distance-dependent in FSO links, and, hence, power allocation schemes that 
take into consideration the distances of the underlying links, can significantly improve system performance. 
Furthermore, it is observed that even the simple sub-optimum power allocations schemes lead to substancial 
performance improvements compared to equal power allocation. Taking into consideration that the parameters 
for these schemes can be easily obtained, based only on the link distances (for most practical FSO applications, 
the link distances are fixed and, thus, are known a-priori at both the transmitter and relays), the proposed 
sub-optimum power allocation can be considered as a less complex alternative to optimum power allocation. 
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VI. Conclusions 

We investigated several transmission protocols for relay- assisted FSO systems without direct link between the 
source and the destination for the Gamma-Gamma channel model. Alternative protocols to the all-active relaying 
scheme were proposed, which activate only a single relay in each transmission slot. Thus, considerable benefits 
in terms of implementation complexity are resulted, since the need for synchronizing the relays' transmissions 
in order for the FSO signals to arrive simultaneously at the destination is avoided. In particular, two different 
types of relay selection protocols were proposed: select-max and DSSC. Select-max relaying offers significant 
performance gains compared to the all-active scheme at the expense of requiring the CSI by all the available 
links. In contrast, DSSC relaying requires less CSI than select-max (only from the links used in the previous 
transmission slot), however it exploits only the two relays with the minimum end-to-end distance. Furthermore, 
based on the derived outage probability expressions, the problem of allocating the power resources to the FSO 
links was addressed, and optimum and sub-optimum solutions that minimize the system's outage probability 
were derived for each considered relaying protocol. Numerical results were provided, which clearly demonstrated 
the improvements in the power efficiency offered by the proposed power allocation schemes. 



Appendix I 

Using the infinite series representation of the Gamma-Gamma pdf ET1 Eqs. (7), (8)] and since Pout,AB 



' AB M fh ( x ) the outage probability for the FSO link between terminals A and B can be obtained 



r h A 

Jo Jh^ 
after some basic algebraic manipulations, as 



P 



outAB ~ r(« A jr(^)^ 



/ 1 

1 KPab+ 1 ) \KbPabPm 



,+i 



\ 



T(l- a AB + f3 AB + 1) 



l\(a AB +l) \h A BpABPM 

T(l + a AB - p AB + 1) 



(52) 



For high values of power margin, i.e., Pm 
approximated by 



oo, the term for I = is dominant and hence (1521) can be 



P 



out,AB 



T Pa~bJ) 



t 



r(a AB )T(P AB ) 



h AB pABP& 



a An$A 



h AB pABPb 



r (1 - a AB + P AB ) /3 AB T(l + a AB - P AB )a, 



(53) 



7 



which can be reduced to (fT9l ), after using the Euler's reflection formula |[T8l Eq. (8.334.3)] and introducing p A 
and q AB . This concludes the proof. 
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Appendix II 

Based on the infinite series representation of the Gamma-Gamma pdf 11211 Eqs. (7), (8)], the pdf of £ 
p h„ n can be written as 



^ W= ' ~*Lw„) ^-^O(^), (54) 

where O (•) represents the least significant terms of an infinite series as x — > oo. By taking the Laplace transform 
of the above equation, the MGF expression of £ m can be obtained as 

a 8 



Hence, the MGF for h SM = P h RmD can be written as 

meS(n) m 



M h ( S ) = s — n A r(a W ) + ° V ^ J ■ (56) 



and, by taking the inverse Laplace transform of (I561 l. an expression for the pdf of h S{n) is obtained as 

<*R m nP Rm n\ qRmD r (lR m n)' r (PR m o-lR m o) 



n 

/„ (x) = m£S(w) ' 7 , x \ meS(n ) -/ + o l xm <rs^ — ) (57) 



K mD p RmD J r(a HmD )r(/3 nmD ) £ q Rm0 -1 / E 9j? 



r( E 9 nmD 

\mSS(n) / 

After some basic algebraic manipulations and keeping only the dominant term, the asymptotic expression in 
(1241) is obtained for the cdf of h s , y This concludes the proof. 
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Fig. 2. Comparison of relaying protocols for a relay-assisted FSO system with d SR . = d R D = 2 km, i G {1, N}. 
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Fig. 3. Comparison of relaying protocols for different relay-assisted FSO configurations: N — 2, d SR = {2, 1.5}, d RD — {1, 2.5} (in 
km) and N = 3, d SR = {2, 1.5, 1}, d RD = {1, 2.5, 3} (in km). 
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Fig. 4. Comparison of power allocation schemes for the relaying protocols under consideration. 
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